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,,^J ' Abstract. We use De Giorgi techniques to prove Holder continuity of weak 

D , solutions to a class of drift-diffusion equations, with L^ initial data and di- 

^^ ■ vergence free drift velocity that lies in Lf°BMOx ■ We apply this result to 

prove global regularity for a family of active scalar equations which includes 
the advection-diffusion equation that has been proposed by Moffatt in the 
context of magnetostrophic turbulence in the Earth's fluid core. 
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1. Introduction 

Active scalar evolution equations have been a topic of considerable study in 
recent years, in part because they arise in many physical models. In particular, 
such equations arc prevalent in fluid dynamics. In this paper we first examine a 
class of drift-diffusion equations for an unknown scalar field 6{t, x), of the form 

dt0+{v-V)e^Ae, (1.1) 

04 ' where v(t, x) is a given divergence free vector field that lies in the function space 

^Ll n L'^BMO-^, i > 0, and a; e W^. In Theorem [Q we prove that weak 
f~^ ■ solutions to (|l.ip are Holder continuous. Note that this result is new for such linear 

O \ parabolic equations with very singular coefficients [ll[18l[19l[25l[27l|3l[35]. We 

^— ^ ■ then use this result to prove in Theorem 13.11 that Leray-Hopf weak solutions of the 

active scalar equation 

dt6 + {u-V)6^ M (1.2) 

^ , divu = (1.3) 

.C^ ; Uj^d^T,,e (1.4) 

are classical solutions. In (|1.4p . the velocity vector u is obtained from 6 via {Tij}, 
a, d X d matrix of Calderon-Zygmund singular integral operators (that is, they are 
bounded L^ H- L^ and L°° ^ BMO) such that did^Tij = 0. Note that in (fTJ]) we 
have used the summation convention on repeated indices, and i, j G {1, . . . , d}. 

Our motivation for addressing the system (|1.2|) - ()1.4p comes from a model pro- 
posed by Moffatt [^^ for magnetostrophic turbulence in the Earth's fluid core. This 
model is derived from the full magnetohydrodynamic equations (MHD) in the con- 
text of a rapidly rotating, density stratified, electrically conducting fluid. After a 
series of approximations relevant to the geodynamo model, a linear relationship is 
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established between the velocity and magnetic vector fields, and the scalar "buoy- 
ancy" 9. The sole remaining nonlinearity in the system occurs in the evolution 
equation for 0, which has the form 

dt6 + {u-y)e = S + KM, (1.5) 

where S* is a source term, and k is the cocfhcient of thermal diffusivity. Here the 
three dimensional velocity u is such that div u = 0, and it is obtained from the 
buoyancy via 

u = M[ei (1.6) 

where M is a nonlocal differential operator of order 1. We describe the precise 
form of the operator M in Section [4l An important feature of this operator is 
the spatial inhomogeneity that occurs due to the underlying mean magnetic field. 
We call (|1.5|) - (|1.6p the magnetogeostrophic equation (MG). We show that the MG 
system satisfies the conditions under which wc prove Theorem l3.1[ and hence obtain 
(cf. Theorem H?T|) global well-posedness for pr3|) - ([LB)) . 

An active scalar equation that has received much attention in the mathematical 
literature following its presentation by Constantin, Majda, and Tabak [5], as a two- 
dimensional toy model for the three-dimensional fluid equations, is the so called 
surface quasi-geostrophic equation (SQG) (see, for example, [2l [6l [7l [TOl [TTl [TBI [29l 
[33] and references therein) . The dissipative form of this equation for which there 
is a physical derivation is 

dt9 + {u-y)9^-{-Ay^^e, (1.7) 

where 

M = V-L(-A)~i/26i = (i?26i,-i?i6i) (1.8) 

and Ri represents the z*'* Riesz transform. It was recently proved by Caffarelli 
and Vasseur [5] that solutions of (|1.7p - (|1.8p with L^ initial data are smooth (see 
also the review article [3]). Well-posedness for (|1.7p - (|1.8p in the case of smooth 
periodic initial data was also obtained by Kiselev, Nazarov, and Volberg [16]. See 
also Constantin and Wu [101 HI] for the super-critically dissipative SQG. 

We note that the magnetogeostrophic equation MG and the critically dissipative 
SQG equation p.7p - (|1.8|) arc both derived from the Navier-Stokes equations in the 
context of a rapidly rotating fluid in a thin shell. For both systems the Coriolis 
force is dominant in the momentum equation. In the case of the SQG equation the 
relation p.8p is derived via a projection of the three-dimensional problem onto the 
two-dimensional horizontal bounding surface. In the case of the MG equation the 
coupling with the magnetic induction equation closes the three-dimensional linear 
system that produces the operators {Tij}, with Uj = diTijO. 

Systems (|1.2p - (|1.4p and (|1.7p - (|1.8p have strong similarities. In particular, they 
have the same relative order of the spatial derivatives between the advection term 
and the diffusive term. Moreover, ii9{t,x) is a solution of (|1.2p - (|1.4p . ihen9\{t,x) = 
9{}?t, \x) is also a solution, and hence L°°{U.'^) is the critical Lebesgue space with 
respect to the natural scaling of the equation. We note that L°° is also the criti- 
cal Lebesgue space for the critically dissipative surface quasi-geostrophic equation 
(|1.7p - (|1.8p . and for the modified surface quasi-geostrophic equation (cf. Constantin, 
Iyer, and Wu [8]). The advantage of system (|1.2p - (|1.4p over the critical SQG equa- 
tion is that the diffusive term is given via a local operator. The tradeoff is that the 
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drift velocity in (|1.2|) - ()1.4p is more singular, i.e., the derivative of a BMO function 
(see Koch and Tataru [17] for the Navier-Stokes equations in BMO^^). 

Our proof of Theorem 12.11 and Theorem 13.11 is along the lines of the proof of 
CafFarelli and Vasseur [2l Theorem 3] for the critical SQG equation. The primary 
technique employed in [2j[TTl[32], and in the present paper, is the De Giorgi itera- 
tion J12j . This consists of first showing that a weak solution is bounded by proving 
that the function max{0 — /i, 0} has zero energy if h is chosen large enough. Then 
a diminishing oscillation result implies smoothness of the solution in a subcritical 
space, namely C", for some a S (0,1). The proof of Holder continuity for solu- 
tions of (|l.ip with V G Lf^BMO^^ does not follow directly either from [2], where 
V G LfBMOx, or from [g, where v G it°°C^-" and a G (0,1). The crucial step 
in the proof of Theorem 12.11 is the local energy and uniform estimates. The main 
obstruction to applying the classical parabolic De Giorgi estimates via an L^'-based 
Caccioppoli inequality (1 < p < cxo), is that v{t^ ■) G BMO^^. In Section[2]wc give 
details as to how we overcome this difficulty. 

Equation (|1.1|) is in the class of parabolic equations in divergence form that have 
been studied extensively, including in the classical papers of Nash [24] , Moser [23] , 
Aronson and Scrrin [1]. Osada [25| allowed for singular coefficients and proved 
Holder continuity of solutions to (|l.ip when v G L'^W~^'°° is divergence free. 
Hence Theorem 12. 11 may be also viewed as an improvement of the results of Osada, 
since if / G BMO{R'^)nL^(R'^), then it does not follow that / G L°°(R'^) (cf. [30]). 
In the same spirit, Zhang [341 I35j and Semenov [2 7) give strong regularity results 
for parabolic equations of the type (|l.ip , where the singular divergence free velocity 
satisfies a certain form boundedness condition. We note that this form boundcdness 
condition docs not cover the case v G L'^BAIO~^, and hence Theorem 12.11 does 
not follow from the results in [27] |34] |35] , and vice- versa. The overall conclusion of 
the body of work on parabolic equations with a singular drift velocity is that the 
divergence free structure of v produces a dramatic gain in regularity of the solution, 
compared to the classical theory (cf. [T5]). 

Organization of the paper. In Section [5] we prove Holder regularity for the 
linear drift-diffusion equation p.ip . with v being a given divergence free vector field 
in the function space L^^ DL'^BMO'^. In Section|3]we apply this result to prove 
that a Leray-Hopf weak solution 6 of the nonlinear active scalar system (|1.2p - (|1.4p 
is Holder smooth for positive time. Since Holder regularity is subcritical for the 
natural scaling of (|1.2p - (|1.4p we can bootstrap to prove higher regularity and hence 
conclude that the solution is a classical solution. In Section 2] we describe an active 
scalar equation that arises as a model for magneto-geostrophic dynamics in the 
Earth's fiuid core. We show that this three dimensional MG equation is an example 
of the general system (|1.2p - (|1.4p . In the Appendix we prove the existence of weak 
solutions to (jl.2p - (|1.4l) evolving from i^(]R'') initial data. 

2. Regularity for a parabolic equation with singular drift 

Consider the evolution of an unknown scalar 9{t, x) given by 

dte + {v ■ V)e = AO (2.1) 

where the velocity vector v{t,x) — {vi{t,x), . . . ,Vd{t,x)) G L^((0,oo) x M'^) is 
given, and {t,x) G [0, oo) x R'*. Additionally let i; satisfy 

d.jVj{t,x)^0 (2.2) 
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in the sense of distributions. We express Vj as 

Vj{t,x) = d,V^j{t,x) (2.3) 

in [0, cxd) X W^, where we have used the summation convention on repeated indices, 
and we denoted Vij = —{—A)~^diVj. The matrix {Vij}f ^i is given, and satisfies 

F„- GL°°((0,oo);L2(M'^))nL2((o,oo);iJi(R'')) (2.4) 

for aU i,j e {f , . . . ,d}. 

Theorem 2.1 (The linear problem). Given 6q G L^(IR'*) and {V^j} satisfying 
(123), let e L°°([0,oo);L2(Kd)) nL2((0,cx));iJi(M'^)) be a global weak solution of 
the initial value problem associated to (|2.ip - (j2.3p . // additionally we have Vij € 
L°°{[to, oo); BMO{M.'^)) for all i,j G {1, . . . ,d} and some to > 0, then there exists 
a > such thate € C"([<o,oo) x R*^). 

In analogy with the constructions in [HIH], the proof of Theorem 12 .11 consists of 
two steps. For <o > fixed, wc first prove that 9 G L°°([to, oo); L°°(E^)). The main 
challenge is to prove the Holder regularity of the solution, which is achieved by using 
the method of De Giorgi iteration (cf. [T2j [T4j [19] ) . Note that for divergence- free 
V G L'^ x; the existence of a weak solution 9 to (|2.ip - (|2.3p . evolving from 9o G L^, 
is known (for instance, see [27] where the more general v G Lj^^ is treated, also [2], 
and references therein). Moreover, this weak solution satisfies the classical energy 
inequality and the level set energy inequalities (|2.6p below. 



Remark 2.2. The conclusion of Theorem 12.11 holds if the Laplacian on the right 
side of (|2.ip is replaced by a generic second-order strongly elliptic operator di{aijdj), 
with bounded measurable coefficients {ay}. 

Remark 2.3. We note that the De Giorgi techniques used here to prove Holder 
regularity for solutions to ()2.ip - p.3p can also be used to prove Holder regularity 
for the problem with a forcing term S on the right side of (|2.ip . In this case we 
consider 5 G L[^ to be an externally given force, with r > 1 + d/2 (cf. [12] )■ 

Remark 2.4. In a very recent preprint, Seregin, Silvestre, Sverak, and Zlatos [55] 
also use De Giorgi techniques to prove Holder regularity of solutions to a parabolic 
equation with drift velocities in L^BMO~^ . 

Notation. In the following we shall use the classical function spaces: L^ 

- Lebesgue spaces, BAIO - functions with bounded mean oscillation, BMO~^ 

- derivatives of BMO functions, H'' - homogeneous Sobolev spaces, and C" - 
Holder spaces. To emphasize the different integrability in space and time we shall 
denote LP([0, oo); L''(M'*)) by LfL« for 1 < p, g < oo, and similarly for LfiJ^ 
and LlBMO^. Also Lf^(/ x S) = L^ {I ; L^ {B)) for any / C M and B C W^. 
The ball in R'' and the parabolic cylinder in R''^^ are classically denoted by 
Bp{xo) = {x G W^:\x- Xf)\ < p} and Qp{to,xo) = [to - p'^,to] x Bp{xo) for 
p > 0. Lastly, we shall write (/ — k)+ = max{/ — k, 0}. 

2.1. Boundedness of the solution. The first step is to show that a weak solution 
is bounded for positive time. 
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Lemma 2.5 (From L^ to L°°). Let 9 e L°°([0, oo); ^^(Rrf)) nL2((0, oo); iJi(M'')) 
he a global weak solution of (|2.ip - ()2.3p evolving from 6q G L^(R'*), where v € 
L2((0, oo); L2(M'^)). T/ien /or all t > we have 

for some sufficiently large positive dimensional constant C . 



Proof. The proof of this lemma is m.utatis- mutandis as in [2l [TT] , and requires only 
the fact that i; is divergence free. The main idea is that since y n- (i/ — /i)+ is 
convex, for all ft. > we have 

dt{e - h)+ - A{9 - h)+ + {v ■ V){0 - h)+ < 0, 

and hence, multiplying by (6 — h)+ integrating by parts, and using that divw = 0, 
we obtain the energy inequality 

*2. 

2, 



mt2,-)-h)+\'dx + 2 \vie-h)+\'dxdt< \{e{ti,-)-h)+\'dx, 

(2.6) 

for all /i > and < ti < t2 < oo. For tg > 0, and i? > to be chosen sufficiently 
large, we define tn ^ t^ — to/2", h^ = H — ii^/2", and 

c„ = sup/ \{e{t,-)-hn)+\^dx + 2 I f \\/{9^h^)\^dxdt, 

where n > 0. The inequality p.6p . the Gagliardo-Nirenberg-Sobolev inequality, 
and Riesz interpolation then imply that 

r ,, < ^ nn(l+4/d) l+2/d 

Letting H = Cc^ jt^ < C||0o||L2(R<i)/to ' ^^r some sufficiently large dimensional 
constant C, implies that c„ -^ exponentially a.sn —> oo, and therefore 0{to,-) < H. 
Applying the same procedure to —6 concludes the proof of the lemma. We refer 
the reader to [5J [TT] for further details. D 

2.2. Local energy and uniform inequalities. In proving the boundedness of 
the solution we only required that v G i? a,, and divw = 0. For the rest of the 
section we use the additional assumption v G L'^BMO"^ . 

Lemma 2.6 (First energy inequality). Let 9 G Lf^L^ n L^H^ be a global weak 
solution of the initial value problem associated to (|2.ip - (|2.3p . Furthermore, assume 
that V,j G L°°((0,cx));BMO(R'^)) for all i,j G {l,...,d}, and ^^ holds. Then 
for any < r < R and ft, G M, we have 

where C ~ C{d, \\Vij\\L°^BMOa:) '■^ ^ fixed positive constant, and we have denoted 
Qp — [to ~ P^yto] X Bp{xo) for p > Q and an arbitrary {tQ,xo) G (0,oo) x M.'^. 
Moreover, estimate (|2.7p also holds with 9 replaced by —0. 



Remark 2.7. Note that from Lemma [2.51 we have that 9 G L^^., and hence the 
right side of (|2.7p is finite. 
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Remark 2.8. The classical local energy inequality (cf. [T4j [TH [25] , see also [2|[TT|) 
does not contain the term \\{d — /i)+||l=° (Qr) on the right, since the velocity field 
V is not as singular as in our case. In this section we prove that since in ()2.7|) the 
exponent 2/(d + 2) of \\{9 — h)+\\[^^^i^Qj^) is "small enough", the De Giorgi program 
may still be carried out to obtain the Holder regularity of weak solutions. 

Proof of LemmalEU Fix /i £ R and let < r < i? be such that to/2 - i?^ > 0. Let 
T]{t,x) e C^((0, oo) X E'') be a smooth cutoff function such that 

< r; < 1 in (0, cx)) x R'', 

r/ = 1 in Qr{xQ, tg), and r/ = in cI{Q'j^(xq, tg) H {{t, x) : t < Iq}}, 

IV77I < j^—, |VVr/| < 19*771 < in QR{xo,to)\Qr{xo,to), 

for some positive dimensional constant C. Define ti = to — i?^ > and let ^2 G [^o — 
r^,to] be arbitrary. Multiply (|2.ip by {9 - h)+7f and then integrate on [^1,^2] x R'' 
to obtain 



dt {{9 - h)l) if dxdt -2 I I djj{e - h)+{9 - h)+Tf dxdt 

+ 11 d,V,j dj {{9 - h)l) if dxdt = 0. (2.. 



The main obstruction to applying classical the de Giorgi estimates (via the L^- 
based Caccioppoh inequality, cf. dHHl]) is that diVij S LfBMO~^, as opposed 
to the case LfW~^°° considered by Osada [25] (see also [27] )■ We overcome 
this difficulty by subtracting from Vij{t, ■) its spatial mean over {t} x Bn, namely 
Vij^BB.{t) (this does not introduce any lower order terms because dxiVij^Bnit) = 0), 
and by appealing to the John-Nirenberg inequality. More precisely, we define 

V^j,R{t, x) = V,j (t, x) - V,,.Bb (t) - V,, {t, X-) - J- / V,, {t, y) dy, (2.9) 

\^R\ J Br 



and note that diVij = diVij.R- Therefore, the third term on the left of (12.81) may 
be replaced by 



•t2l 




^^V^J,R dj {{9 - h)l) rf dxdt. 
1 



We integrate by parts in t the first term on the left of (|2.8p . and use 77(^1, •) = 0. 
The second term we integrate twice by parts in Xj , and the third term on the left of 
we integrate by parts first in Xj (and use dj{diVij^R) = di{djVij) = djVj ~ 0) 
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and then integrate by parts in Xi, to obtain 



i / ieihr)^hjir,ih,-rdx + 



rf \\7{e-h)+fri^ dxdt 



[e - h)\ridtri dxdt + / (6- h)ldj{rjdjj]) dxdt 
I 7 %.r{0 - h)ld,iTjdjTj) dxdt 

Vij,R^^iO - h)+{e - h)+Tjdjr] dxdt. (2.10) 

Using the bounds on the time and space derivatives of 77, the fact that 77 = 1 on 
Qr, t2 < to, the Holder and e- Young inequahties, we obtain from (|2.10p 



/ {e{t2,-)-h)ldx + 2 f 'f \V{e-h)+\^Tf dxdt 



(R^ lif " '^^ ''"' + (^ J J,. ' ^- 



< 77r—v2 II (^ - '^)l dxdt + Tj^—K2 1 1 l^y.fll (^ - '^)+ dxdt 

C 



, , \V{e^h)+\^Tfdxdt+ ^ II \V,j,R{'{d^h)l dxdt. (2.11) 

tiJR'' [^^^) JJQr 

After absorbing the third term on the right of (|2.11l) into the left side, we take the 
suprcmum over ^2 G [^o ~ ''^ > ^0] : to obtain 

- (i?^"^^ ^^^WllAQ.) + (]^ // l%«l (^ - '^)+ dxdt 



C 
- I 

'Qr 



^^_^)2 // \V^,,R\'{e-h)idxdt. (2.12) 



As a corollary of the celebrated John-Nirenberg inequality (cf. [TH [3D]) we have 
that for any fixed R > 0, t & [to — E? , to], and 1 <p < 00, 

\\%At.-)\\LnBn) = \\V.,{t,-)^V^j,BAt)\\L.iBR) 
<C\\V,,{t,-)\\BMO{R-)\BR\"'', 

where C = C{d,p) > is a fixed constant (recall that C{d,p) — > 00 as p — !• 00). 
The fact that V^j e L°°([io/2,oo); BMO(R'')) implies that for all i G [to - i?^to] 
we have 

m,M{tr)\\LHBn)<Co\BR\'/P (2.13) 

for a positive constant Co = CodlK j ||L°°([to/2,oo);BMO(R'i)), ti.p). We fix < £ < 2 
to be chosen later, and using (|2.13p and the Holder inequality we bound 

rto 
\2 



|Kj,fl,| (0 - h)\_ dxdt = / / \V,jAt,x)\ (6 - hy_^{t,x) dx ) d< 



<Co\Br\''' / ||(0(t,.)-/i)+|li4/(.-.,(B«)di. 

Jto-R? 
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Using the interpolation inequality \\f\\Lp < C'||/||/2''||/||i=o , with p = 4/(2 — e), 
we obtain from the above estimate that 

V^y,i?J [0 - h)l dxdt 

Qr 

to 

dt 



<Coi?^(^+^)/^||(0-/^)+||^r(Q,)||(^-/.)+|li^^(Q,). (2.14) 

Similarly, from (j2.13p . the Holder inequality and L^ interpolation, we obtain 

\V.j,R\\d - h)l dxdt < CoE^('^+2)/2||(0 - /^)+llir (Q,)ll(^ - h)+\\l^jQ^y 

(2.15) 

Combining estimates (|2.12p with p.l4p . (|2.15p . and the Holder inequality, we con- 
clude that 

fi ne(d+2)/2 

The proof of the lemma is concluded by letting e = 2/((i + 2) in (|2.16p above. D 

By applying the Holder inequality to the right side of (j2.7p we then obtain: 
Corollary 2.9. Let 6 he as in Lemma \2.6[ Then we have 

r<f?d+2 

m /^)+lli~i.(Q.) + l|V(0 - M + ili?,^(Q.) < (]^T7)^II(^ - hMl-JQ.y (2-17) 
for some positive constant C ~ C{d, ||Vij ||l°°smOx)- 

We now fix a point {to, xq) G (0, oo) x W^ and we prove the Holder continuity of 
9 at this point. Throughout the following we denote by Qp the cylinder Qp(tQ,xo), 
for any p > 0. 

The following lemma gives an estimate on the supremum of on a half cylinder, 
in terms of the supremum on the full cylinder. A similar statement may be proven 
for -0. 



Lemma 2.10. Let 9 be as in Lemma \2.6\ Assume that Hq < supg 6, where rg > 
is arbitrary. We have 

for some positive constant C ~ C{d, ||T^y ||l°°sa/Ox)- 

The above estimate differs from the classical one cf. [TH Theorem 6.50] in that 
the power of \{9 > /io}nQ,ol/IQr-ol is 1/(2^ + 4) instead of 1/(^ + 2). However, the 
key feature of (|2.18p is that the coefficient of (supg^ 9 — ho) does not scale with 
tq. It is convenient to introduce the following notation: 

• A{h,r) = {9> h}nQr 

• a{h,r) = \A{h,r)\ 



. bih,r) = \\{9-h)+r, 
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• AI{r) ~ SUpg^ 9 

• m{r) = infg,, 9 

• osc(Q) = supg 9 — inf Q 9 

Proof of Lemma [2Jd[ Let < r < i? and < h < H. We have 

b{h, r) = \\9- /^lli. j^(,,,)) > 11^ - /^lli? jA(i/..)) > (^ - hfa{H, r). (2.19) 

Let r]{t,x) G C^{R x R'') be a smooth cutoff such that 77 = 1 on Qr, 77 = on 
Q'(r+Fi.)/2 ^ {^ — ^0}, and IV77I < C/{R — r) for some universal constant C > 0. 
Then, by Holder's inequality and the choice of 77 we obtain 



b{h,r) = \\{9 - /^)+|li.^(Q,) < a(/7,r)2/('^+2)||(0 - /i)+||^,,,+,„. 



{Qr- 



<a{h,rf/^^+^^\\7j{9-h) "2 



+ l'L?<,^+">''''((-oo,to)xR^)- 

(2.20) 
Using the Gagliardo-Nircnberg-Sobolev inequality and Riesz interpolation 

||/|lL2(d+2)/d((_oo,to)xR'i) - ^\\f\\L^Ll{{-oo,to)xR'') + '^ll^/llL?_^((-oo,io) xR'i) ' 

estimate (|2.20p implies that 

6(/7,r)<Ca(/7,r)2/('^+2)(h(0-/7)+||i^^.((_^^,^),^.) 

for some positive dimensional constant C. Using the first energy inequality, i.e., 
Lemma l2.6i and the Holder inequality, we bound the far right side of the above and 
obtain 

6(/.,o<Ca(;7,.)^/(^-)^^ll(.-/7),||-X)')||(.-/7),y^ 

< Ca{h, rf/^^^^) jR^bih, i?) 1- V(<i+2) II (^ _ /^)^ ||2/(.+2) ^ (2 21) 

for some sufficiently large positive constant C = C{d, \\Vij\\L°^BMO^)- By combining 
estimates p.l9p and p.2ip we obtain the main consequence of Lemma 12.61 that is 

^(^'-) ^ (ff_fe)W)(^_,). Hh,Rr'^^'-'^iO-HU\\l/J,^^l^. (2.22) 

The above estimates give the proof of the lemma as follows. Let r„ = ro/2 + 
ro/2"+^ \ ro/2, hn = hoc - {hoo - /io)/2" /■ /loo, and 6„ = 6(/i„, r„+i), for all 
71 > 0, where t'o and h^ are as in the statement of the lemma, while /iqo > is to 
be chosen later. By letting H = hn+i, h = hn,r = r„+2, and R = r„+i in (|2.22p . 
we obtain 

C(A/(ro) - /lo)^/('^+^) r,n(2+4/(rf+2)), l + l/(d+2) (r, r,o\ 
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by using 

||(6'-/i„+i)+||ioo(Q )= sup e~K+i 

< sup 9 — hn+i < svLp0 — ho = M(ro) — ho 

which holds since M^r^i) < Af (ro) and h„ > hg. Let B = 2"*+2(''+2). We choose h^o 
large enough so that 

then by induction we obtain from p.23p that &„ < bo/B^\ and therefore 6„ — > as 
n -^ oo. This implies that supg 9 < hoo- A simple calculation shows that if we 
let 

ftoo - /lO + (^1^5)74 (.^-^Sj 

then (|2.24p holds. Lastly, bo — 6(/io,3ro/4) may be bounded via (|2.2ip and the 
Holder inequality as 

bo < Caiho,ro)'/^^+^^4+^iMiro) - hof . (2.26) 

The proof of the lemma is concluded by combining supg < ^oo with (|2.25p 
and (|2.26p . From the above proof it follows that inequality ()2.18p also holds with 
9 is replaced by —6*. D 

As opposed to the elliptic case, in the parabolic theory we need an additional 
energy inequality to control the possible growth of level sets of the solution. 

Lemma 2.11 (Second energy inequality). LeA 9 G LfV^^ n L^H]^ be a global 
weak solution of the initial value problem associated to (|2.1[) - (j2.3p . Furthermore, 
assume that Vij € Lf°BMOx for all i,j e {!,..., d}, and p.4|) holds. Fix an 
arbitrary xo S K."^, let h gM., < r < R, and < ti < <2- Then we have 

11(^(^2, •)-/^)+lli.(B„) 

< \mt„ .) - /.)+||i.(^,) + ^^^^l^\\{9 - M+lli^j(*,,.)xB.) (2.27) 

for some sufficiently large positive constant C — C(d, ||14j [|L°°BAfOx)7 where we 
have denoted Bp — Bp{xo) for p > 0. 



Proof. Note that by Lemma [2.51 we have that 9 g L^ and hence the right side of 
^ftjn\ is finite. Let 77 G C^iW^) be a smooth cutoff such that 77 = 1 on B^, r; = 
on Bfj,, and |Vr/(a-)| < C/{R-r), for ah x G R'^, for some constant C > 0. Multiply 
p.ip by ri^{9 — h)+ and integrate from ti to ^2 to obtain 



dt(i9 - h)+) jfdxdt -2/7 djj{9 - h)+{9 ~ h)+ri^dxdt 

i{i9-h)+)\ 



*2/- / ^2 

tl 

2 f 

diVij^R dj ( {9 - h)+ ] ri^dxdt, 



d,V,j dj[{9-h)+\ n^dxdt 



*2 /• , ^ 2 

2, 
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where, as in p.9p . we have denoted Vij^R{t,x) = Vij{t,x) — t^— r /g Vij{t,y) dy. 
After integrating by parts we get 

^2 „2 . , / / |V7//) h\ |2^2, 



{eit2,-)-h)lv dx+ / \S7ie-h)+\''r]''dxdt 

< f {e{h,-)-h)l7fdx 

+ C fj {d- h)l (\dj{'nd,r])\ + \V,j.R\mvdjv)\ + \V^j.R\^\d,vd,v\) dxdt. 

We bound the right side of the above estimate as in (j2.14p and (|2.15p to obtain that 

||(0(i2,•)-/^)+||i.(B.) 



< ||(0(ti, •) - /i)+||i.(s,) + {^i^_,y W h)AY{j{W - h)A\ 



L? 



Letting e = 2 in the above estimates concludes the proof of the lemma. The 
corresponding statement for —9 also holds. D 

The use of the second energy inequality is to bound |{0(t2, •) > H} D Br\/\Br\, 
whenever |{6'(ti, ■) > h} n Br\/\Br\ < 1/2. More precisely, we have the following 
lemma. 

Lemma 2.12. Fix ko = (4/5)^/'', let no >2 be the least integer so that 2"o/(2"'' - 
2) < y/6/5, and let (5o = (1 — ko)^/(12CoKq), where Cq is the constant from (j2.27p . 
Forti,R> 0, if 

\{0{ti,-)>h}nBr\<^\Br.l (2.28) 

then for all ^2 G [^i, ^i + Sor"] we have 

mt2r)>H}nBn\<^-\BRl (2.29) 

where r ^ kqR, M == sup(t^ t^+^^^s)^^^ 6*, m = val(t^,t^+8aB?)^BRQ , h = {M + 
m)/2, and H ^ M - {M - m)/2"o. 

Proof. For i2 G [ti,ti+6or'^], we obtain from the second energy inequality (cf. (|2.27p ) 
that 

||(0(t2, •) - /i)+lii.(B,o ^ 11(^(^1' ■) - f')+\\hiB.) + ^"fj^'y^^ lK^ - f'^+WlrAQ.) 

< mt,, .) - /i)+||i.(5,) + J^-^{M - h)\ (2.30) 

where Q2 ~ {ti,ti + 6qR'^) x Br. The left side of the above estimate is bounded 
from below as 

||(0(i2,-)-/i)+lli2(B,) > ||(0(i2,-)-/i)+lli.(B.nWt„.)>//}) 

>{H~hf\{e{t2,-)>H}nBr\. (2.31) 

From (j2.30p . p.Sip . and the Holder inequality, we obtain after dividing by \Br\ that 

mt2,-)>H}nBr\ ^{M-h)^R'' f\{e{ti,-)>h}nBR\ , CoV' 



Br\ - [H-hfrd \ \Br\ [l-r/RYR^ 
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Noting that by construction (M - h)/{H - h) ^ 2"n/(2"° - 2) < ^/&Jh, and 
recalling that r/R ~ kq ~ {4/5)^/'^, wc obtain from the previous estimate and the 
assumption of the lemma that 

mt2,-)>H}nBr\ ^3 f\{e{h,-)>h}nBR\ , CoSoi^l 



\Br\ -2\ \Bn\ (1-k; 

concluding the proof of the lemma. D 

2.3. Holder regularity of the solution. Wc now have all necessary ingredients 
to conclude the Dc Giorgi argument for proving Holder regularity of the weak 
solution. 

Recall that since divu = 0, by Lemma [231 we have that 6 e L°°i[to, oo); L°°(R'^)) 
for any to > 0. Moreover, if V^-j G L°°{[tQ,oo); BMO{R'^)) for some to > 0, we 
obtain the energy inequalities of Lemmas 12.61 and 12.111 In turn, these inequalities 
give control for the growth of the supremum on doubling cylinders (cf. Lemma [2.10p , 
and for the growth of level sets of the solution (cf. Lemma I2.12p . The rest of the 
proof follows as in ^^, but we give a sketch for the sake of completeness. 

Proof of Theorem \2.1\ The proof of the theorem is based on showing that there 
exists 7 £ (0, 1) such that osc((5i) < 7 osc{Q2)- The key observation is that if 7 is 
independent of i?, this estimate implies the Holder regularity of the solution, where 
the Holder exponent a G (0, 1) may be calculated explicitly from 7. 

Fix kq, ^Oi '^Oi ^1 TTi-, h, H, r, and R as in Lemma 12.121 for the rest of this proof. 
We also fix two cylinders Qi ~ [ti,ti + Sor'^] x Br, and Q2 = [ti,ii + 6oR^] x Bn, 
where we recall that ti > Q and R > are arbitrary. 

Recall that h = (infgj 6' + supQ2 6')/2. Without loss of generality we may assume 
|{6l(ti, ■)>h}n Br\ < \Br\/2. Otherwise , letting h' = (infQ,(-6') + supQ^(-6'))/2 
we have |{-6'(ti, •) > h'} n Br\ = |{6'(ti, ■) < h} n Br\ < |S,.|/2, and we work with 
—9 instead of 0. 

For 71 > no, we define i/„ = AI — [M — m)/2", and note that H ~ Hn^ < Hn /" 
M . We also let w be 6* truncated between levels i?,i-i and iJ,i, namely 

0, & < Hn-i 

w = min{6i, Hn} - min{6i, H„-i} = { 9- Hn-i, H^-i < 8 < H^ 

Hn — Hn-1, Hn < 0. 

Since \{9(ti, ■) > h} n Br\ < \Br-\/2, by Lemma [^1^ for every t e [ti,ti + Sor'^] we 
have 

\{wit, •) = 0} n Br\ = \{e{t, •) < Hn-i} n Br\ > \{e{t, ■)<H}n Br\ > ^\Br\. 

By the above estimate and the Poincare inequality we obtain 

\w{t,-)\dx <Cr \Vw{t,-)\dx 

I Br- J Br 

for all i e [ii , ti+Sor"^] , where C = C((i) is a universal positive constant. Integrating 
the above estimate in time over [ii,ti + ^o''^] and using the Holder inequality we 
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get 



// \w\dxdt<Cr \Vw\dxdt 



< Cr|{H„_i < < i/„} n Qi\^/^\\V{e - i/„_i)+L2^(Q^). (2.33) 
We bound the far right side of (|2.33p by using Corollary 12.91 to obtain 

\w\dxdt < Cr|{iJ„_i <e<H^}n Qi\'/^\\V{0 - H^-i)+\\Lr (Q2)^i^ 
Qi ■'^ H-r 

< C-^^^\{Hn-i < e < Hn} n Qi|i/2|Q2r/'(A/ - Hn-i) (2.34) 

1 - KO 

The left side of (|2.34p is bounded from below as 



w\dxdt> \w\dxdt>{H,,-H„_i)\{0>H„}r]Qi\. (2.35) 

Qi J JQin{e>H„} 

By combining and squaring estimates p.34p and (|2.35p we obtain 

< C\Q2\ {\{e > Hn-i} r\Qi\-\{6> Hn} n Qil) , (2.36) 
where we used the faet that, by construetion, {M — Hn-i) / {Hn — Hn-i) = 2. Hence, 
^ \{o > Hn} n Oip < C\Q2\\{0 >Hn„}nQi\, 

n>no4-l 

and since the sequence |{^ > Hn} CiQil is decreasing, we obtain 

C\Q2\'^'\{0>H}nQi\^/^ 



\{0>Hn}nQi\< 



(n-no)i/2 



for all 71 > 710 + 1. By Lemma 1^1^ we have that \{e > H}r\Qi\ < 7\Qi\/8, and 
therefore the above estimate implies 

\{e>Hn}nQ,\< ^'' (2.37) 

(71 -no)^/-^ 

where we have used that r = kqR, and kq = Ko{d). By Lemma 12.101 the fact that 
So < 1, and the estimate (|2.37p we obtain 

Qi V '^ / 

^^"+ (71-7JV(4.+B) (^^-^")' 

for some positive constant C = C{d, \\Vij\\L^BMOo:)i which is independent of r. 
Therefore there exists a sufficiently large rti = ni{d, WVijUl'^bmo^) > ?^o + 1 such 
that 

SVip0<Hn, +-{M -Hn,). 
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Recalling the definition of i/„, m, and M, a simple calculation shows that the above 
estimate implies 

osc(Qi) = sup0 - inf < H„, - m + 1{M -H„,)^(l- — ^ ) (Af - m) 

= (l - ^) («UP^ - i^f ^) = ^ °«c(Q2), (2.38) 

where 7 = 1- l/2"i+2 e (0, 1) is independent of r. Recall that in ([OS)) we 
have Qi = [ti,ti + SquIR"^] x Bko-rX-^o) and Q2 = [^1,^1 + SqR^] x Br{xo), with 
Ko7'^o fixed positive constants, and i? > arbitrary. This classically implies Holder 
continuity of 9 at the arbitrary point {ti.xo) £ (0,oo) x W^, concluding the proof 
of the theorem. D 

3. Global regularity for a nonlinear parabolic equation 
We address the global regularity of solutions to the initial value problem 

dtO- M+{u-y)e = Q (3.1) 

divu = (3.2) 

uj = d.,T,,9 (3.3) 

e(0,-) = ^o, (3.4) 

where {Tij}f ^^i is a matrix of Calderon-Zygmund singular integral operators such 
that didjTijf ~ for any Schwartz function /. As an elementary example, if d = 2 
we may consider Tn = T22 = 0, and T12 = — T21 = T, for some Calderon-Zygmund 
operator T (for instance T = Ri, a Riesz-Transform) . In this case the velocity 
would be M = V-^T9. When rf = 3, a physical example of such a matrix {Tij} arises 
in the MG system (cf. Section |4] below) . 

Theorem 3.1 (The nonlinear problem). Let 9q € L^(R'') be given. A Leray- 
Hopf weak solution 9 e L°°([0, 00); L2(R^)) n L^ {{0, 00); H^{W^)) of l^ -UKM. 
evolving from 9o, is a classical solution, that is 9 ^ C°°((0, 00) x W^). 

Lemma 3.2 (Boundedness). A Leray-Hopf weak solution 9 of (|3.1[) - p.4[) is 
hounded fort > 0, i.e., 9 G i°°([to, 00); L°°(M'^)) for any to > 0. 

Proof. The proof of this lemma is the same as the proof of Lemma [^31 (cf. pilllj). 
and only uses the fact that divw = 0, where u 6 ij ^((0, 00) x M''). D 

Since 9 G i^, it follows from the Calderon-Zygmund theory of singular integrals 
that T^j9 ~. V^j G £°°([<o, 00); BMO{W^)), for any io > 0, where i, j G {1, . . . , d}. 
Therefore, we may treat (|3.ip as a linear evolution equation (see also [H [TTj). 
where the divergence-free velocity field u is given, and u G i^((0, 00); L^(M'')) n 
i°°([io,oo);BMO-i(R'^)), for any to > 0. This is precisely the setting of Theo- 
rem [5TT] for the linear evolution equation. Hence Theorem 12 . II can be applied to the 
nonlinear problem to give Holder regularity of the solution. Therefore we obtain: 

Lemma 3.3 (Holder regularity). A Leray-Hopf weak solution 9 of (|3.ip - p.4p 
is Holder smooth for positive time, i.e., for any to > 0, there exists a > such that 

6'GC"([io,oo) xR'*). 

Lastly, since the Holder regularity is sub-critical for the natural scaling of p.ip - 
p.4p one may bootstrap to prove that the solution is in a higher Holder class: 
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Lemma 3.4 (Higher regularity). Let 9 e L°°([to, oo); C"(R'')) be a Leray-Hopf 
weak solution of the initial value problem associated to (j2.1|) ~ (|2.3p . with a G (0, 1). 
Then 9 e L°°([ii, oo); C'^+^{R'^)), for any ti > to, for some S G (0, 1). 

For 1/2 < a < 1, the proof is the same as the proof of higher regularity for the 
modified surface quasi-geostrophic equation [H Theorem 2.2] (see also ^0\ Theorem 
3.1] for the supercritical quasi-geostrophic equation). These elegant proofs use the 
natural characterization of Holder spaces in terms of Besov spaces, and energy 
inequalities at the level of frequency shells. 

For < a < 1/2, the C" smoothness of 9 is weak relative to the roughness of the 
velocity u, and it is therefore necessary to modify the techniques of [TOl [8] for the 
proof of higher regularity. In [13] we give the details of this modification which uses 
the extra information that u G Lf ^ and employs estimates in the Chemin-Lerner 
(cf. [5]) space-time Besov spaces. 

We give a very brief outline of the proof of lemma 3.4 in the two ranges for alpha 
and refer the reader to [TOl |8] and [13] for detailed estimates. 

Proof of Lemma \3^ Let Bt, „ be the classical homogenous Besov space (cf. [SlfTO]). 
and recall that L°° n B;^ ^^ = C* is the Holder space with index s. The proof of 
the lemma in the case a S (1/2, 1) is based on first noting that if 9 is as in the 
statement of the lemma, then 9 £ i°°([io, oo); Bp^oo), where a^, = (1 — 2/p)a, and 
p e [2, cx)) is fixed, to be chosen later. Then, for j € Z fixed, we have 

~\\^A\l, + J \A,9r''A,9{-A)A,9 = -J \A,9r^A,9A,{u ■ V9). (3.5) 

Upon integration by parts (see also [J), the dissipative term is bounded from below 

j \A,9r'A,9{~A)A,9 dx > ^^—^\\A,9\\1,, (3.6) 

where C'(d,p) > is a constant depending on the dimension and p. The main 
difficulty lies in estimating the convection term. This is achieved in [8j [10] by using 
the Bony paraproduct formula, the Holder inequality, the Bernstein inequalities, a 
commutator estimate, and the fact that ||u||(;^op-i < C||6'||(7"p. The latter holds 
since Uj = diTij9 and the fact that Calderon-Zygmund operators are bounded on 
Holder spaces. If a,, < 2 these operations give 



/ \Aj9Y'-'^Aj9Aj{u ■ \79) dx 



<^2(2-2".b||0||^„^||0|| (3.7) 



Combining p.5p - p.7p . using the Gronwall inequality, and then taking the supre- 
mum in j gives that 9 G i°°([ti, oo); i3p^(R'^)) for any ti > to- Using the Besov 
embedding theorem we obtain that 9 G L°°([ii, oo); -600,00" (R'^)), for any ti > to, 
where €p = (4q; + d)/p < (4-1- d)/p. Letting p > {4 + d)/{2a — 1) concludes the 
proof of the lemma in the case a G (1/2, 1). 

In the case a G (0,1/2] the proof is based on proving that the additional infor- 
mation 9 G L^([<i,<2];-ff"^), implies 9 G L^([<i,i2]; ^^^^ '^) for some large enough 
p > 2, and for any ^2 > ii. This is achieved by using the smoothing effect of the 
Laplacian on high frequencies of 9, so that we need to work in the space-time Besov 
spaces introduced by Chemin and Lerner (cf. [5]). By the endpoint Sobolev embed- 
ding theorem we thus obtain that W9 G L^([ti,t2]; ^So 1) C I/^([ii,t2];-^°°)- From 
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here, standard energy estimates imply that 9 E L°°{[t[,t2];H^"') for all m > 2, 
and t'l G [ti , ^2] , concluding the proof of the lemma after applying the Sobolev 
embedding H"' C C^'^ with m > 1 + d/2. We refer to [T^ for details. 

D 



Proof of Theorem \3.1[ The existence of a global in time Leray-Hopf weak solution 
of ()3.ip - (|3.4p . evolving from 9o G L^, is proven in Appendix [Al The argument is 
to construct solutions to an approximate system, and then to pass to the limit in 
the weak formulation of the problem, using the Aubin-Lions compactness lemma 

(cf. m)- 

The proof of Theorem 13.11 now follows from Lemmas 13.21 13.31 13.41 For any 
/3 G (0,1), after finitely many applications of Lemma 13.41 the solution is shown 
to be in L°°{[to,oo);C^'^^{M.'^)), for any to > 0, and is hence a classical solution. 
Higher regularity is standard. D 

4. Global regularity of the MG system 

There is a vast literature studying mathematical models for the Earth's dynamo 
(see, for example Glatzmaicr, Ogden, and Clune |15j and references therein). How- 
ever, at present, no computational dynamo model can encompass the fine scale res- 
olution required to simulate the turbulent processes believed to exist in the Earth's 
core. It is therefore reasonable to examine models that are simpler than the full sys- 
tem of PDE governing rotating, convective, magneto-hydrodynamic flows, but that 
retain some of the essential features relevant to the physics of the Earth's core. One 
such model for magnetostrophic turbulence was recently proposed by Moffatt [22] . 
He postulates that the magnetic field B{t,x) in the core consists of a mean part 
Bq, which results from dynamo action and can be considered as locally uniform 
and steady, and a perturbation field 6(t, x) induced by the flow w(t, x) across -Bo- 
It is assumed that the scale L of convective turbulence lies in the range y/fi ^ 
L <C Tj/V, where V is the average magnitude of the upward buoyant velocity, 51 
is the angular velocity of the Earth, and 77 is the magnetic diffusivity of the fluid 
medium. This assumption implies that the Rossby number V/flL and the magnetic 
Reynolds number VL/rj are both small. The turbulent Reynolds number in the core 
is expected to be very large. The dominant terms in the three dimensional equations 
of motion and the induction equation give the following linear system 

21763 X li = -VP + (Bo • V)6 - % (4.1) 

= (Bo ■ V)w + r/A6 (4.2) 

divM^O (4.3) 

divfe = 0, (4.4) 

where P{t, x) is the sum of the fiuid and magnetic pressures, 9{t, x) is the buoyancy 
field (e.g. perturbation of the temperature), and g is the gravitational acceleration. 
We use Cartesian coordinates in the reference frame rotating about the axis 63 = 
(0,0,1). 

Equations (|4.ip - (j4.4p establish a linear relation between the variables u{t,x), 
h{t,x), and 9{t,x). The sole remaining nonlincarity from the full convective MHD 
system occurs in the advection-diffusion equation for the buoyancy 9{t,x): 

dt9 + {u-\7)9 = S + kA9, (4.5) 
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where S* is a source term. The difFusivity k in the core is very small, hence the 
nonlinear advection term is dominant and cannot be neglected. 

The system (|4.ip - (|4.5p gives an active scalar model for magneto-geostrophic dy- 
namics, which we call the MG equations. As Moffatt observes, (|4.ip - (|4.5p has some 
similarities with the dissipative Burgers equation, but it has a clearer physical basis 
and the velocity uit, x) is three-dimensional. We remark that the system has closer 
similarities to the surface quasi-geostrophic equation (SQG), which is also derived 
in the context of a rapidly rotating system dominated by Coriolis' force. However, 
the operator that connects u and via (|4.ip - (|4.4p has features that are distinct 
from the analogous operator in the SQG system as we shall now discuss. 

For simplicity we will examine (|4.ip - ()4.4p in the case where Bq is a vector that 
is constant in magnitude and direction in the plane perpendicular to 63. We write 

We assume that gravity acts parallel to the axis of rotation, i.e. g = e^. With 
these assumptions we are examining a local tangent plane model for the Earth's 
fluid core that ignores the sphericity, but retains the essence of the mathematical 
structure of the active scalar equation (|4.5p . with u constructed from 9 via (|4.ip - 
(|4.4p . Manipulation of the linear system (j4.ip - (|4.4p gives, in component form, 

ui = D-^ {-2nd2P - TdiP) (4.6) 

U2 = D-^ {2VLdiP - Td2P) (4.7) 

(4.8) 

P, (4.9) 







d3U3 = D~^TAhP 






d3e^{T^AHD~'+d33) 


where the operators F, 


D, 


and Ah 


arc defined as 






F== 


_^(_A)-ia22 
v 






D = 


An^ + T^, 






Ah 


== 9ii +922, 



(4.10) 

(4.11) 
(4.12) 

where x ~ {xi,X2,X3) e M^ x T. We note that a more general choice of the mean, 
steady, locally uniform magnetic field Bq or of the gravitational vector g results in 
the same structure of the leading order terms. It is the anisotropy that is produced 
by Bq that is a distinctive and crucial feature of the MG system. 

The operator D given by (|4.1ip is invertible since its Fourier symbol does not 
vanish on M^ x Z, justifying the use of _D~^ . In order to uniquely determine M3 and 9 
from (|4.8p and ()4.9p . we restrict the system to the function spaces where 9 and U3 are 
periodic in the xa-variable, with zero vertical mean, i.e. L 9 dx^ = L U3 dx^ = 0. 
In fact, without such a restriction the system is not well defined. We integrate (|4.9p 
and use the zero-mean assumption to obtain 

9 = A[P]. (4.13) 

where A is formally defined as the Fourier multiplier with symbol 



4»^fc||fcp-H(/3V??)^fc.4 



A{k,,k2M - .. ..o^Ju ,V;2/".2^ (4-14) 



for all fcs f^ (by our vertical mean- free assumption), where k = (fci,fc2,fc3) S 
M? X Z. Therefore A is invertible on the space of functions with null a;3-average. 
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Note that d^AlP] = {T^AhD-^ + d33)P in the physical space. We now use (gj])- 
()4.8p to represent ui,U2, and u^ in terms of 9: 

ui ^ D-\-2nd2 - rdi){A-^[9]) = Mi[9] (4.15) 

U2 = D-\2ndi - rd2){A-^[9]) = M2[9] (4.16) 

U3 = {D-'TAhKD-'TAh + d3:i)-'[9] = M3[9]. (4.17) 

To investigate the properties of the operator M = (Mi, M2, M3), we note that it is 
a vector of Fourier multiphers, with exphcit Fourier symbols given by 

for all A;3 7^ 0. Since by assumption 9{ki, fc2, 0) — u{ki, fc2, 0) = 0, in order to have 
a uniquely defined symbol M{k) on all of M^ x Z, without loss of generality we may 
let Mi(fci,fc2,0) = M2(fci,fc2,0) = 0, and M3(fci,fc2,0) = M3(fci,fc2,l). Note that 
Uj ~ Mj [9] is defined via the inverse Fourier transform from 

Uj{k) := Mj{k)9{k), for all fc G M^ x Z, (4.21) 

for all j G {1, 2, 3}. Also, since divu = 0, we have that k ■ M{k) = 0. 

When the frequency vector k = {ki,k2,k3) has components such that fci < 
max{fc2, fc3}, then the symbols Mj arc bounded for all j G {1, 2, 3}. However this is 
not the case for "curved" regions of frequency space where ^3 — 0(1), fc2 — 0{\ki\'^), 
where < cr < 1/2, and |fci| ;» 1. In such regions the symbols (|4.18p - ()4.20p are 
unbounded, since as |A:i| -^^ 00 we have 



iA./i(fci,|fcir,i)i«|fcir, iM2(fci,|fcir,i)i«|fcii, |A/3(fci,ifcir,i)i«ifci 



\2a 



where a G (0, 1/2], and we write a ?» fe if there exists a constant C > such that 
a/C <b <Ca. It follows from gUf-gSQl) that 

|Mj(/c)| <C*|fc| (4.22) 

for aU k e M."^ xZ, and all j G {1,2,3}, where C* = C*(/3,ry,f)) > is a fixed 
constant. From the previous remark it is clear that along certain curves in frequency 
space the bound (|4.22p is sharp. 

Wc now prove that the active scalar equation (cf. (|4.ip - (|4.5p with S" = 0) 

dt9 + {u ■ \7)9 = kA9 (4.23) 

divu = (4.24) 

u = M[9] (4.25) 

with M given by (|4.15l) - (|4.17p . or cquivalently by its Fourier symbol (|4.18p - (|4.20p . 
satisfies the conditions of the abstract problem studied in Section [3] First note that 
we can write 

Uj = M, [9] = 9,Ty- [9] = d.,V,j , (4.26) 
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where we have denoted 

Ty = -d,i-A)-Hlj. (4.27) 

By (|4.22p we have that \Tij{k)\ < C* for all fc £ R^ x Z, and hence it follows directly 
from Planchcrel's theorem that T^ : L'^{R'^ x T) ^ L'^{R'^ x T) is a bounded map. 
It remains to prove that T,j: L°°{R'^ x T) i-^ BMO{M.'^ x T) boundedly. This 
reduces to proving that Nj = (— A)~^/^Mj : L°° H> BMO is a bounded map, 
since Riesz-transforms arc bounded on BMO. The later holds because Nj is a 
pseudo-differential operator of order (cf. [211 EH EO]). The main idea is that 
one may extend Nj from M^ x Z to a symbol iV' defined on R? such that they 
agree on M^ x Z, and such that TV' is the symbol of a classical Hormandcr-class 
pseudo-differential operator of order (cf. Stein [30]). More precisely, let N'Ak) = 
Alj{k)/\k\ for all fc £ R'^ with |fc3| > 1, while for \k^\ < 1, replace the denominator 



4f72fc2|^|2 _^ (/32/^)2^4 by ^j^^ quantity ifl^ipik^fikl + kj + ip{k-if) + {P'^/rifk: 
the definitions (|4.18p - (|4.20p of Mj(fc). Here Lp{-) is C°° smooth monotone increasing 
function that coincides with the identity on jfeal > 1/2, and is constantly equal to 
1/2 on Ifcsl < 1/4. This construction ensures the smoothness of the symbol near the 
origin, while the bound \d'^iS['Ak)\ < Cq(1 -|- |fc|)~l"l follows by inspection. To close 
the argument, note that the operators Nj and N', differ by a compact operator in 
the symbol class S'~°"(cf. [21][26| and references therein). This concludes the proof 
of the boundedness of T,j : L°° H- BMO. 



The abstract Theorem l3.1l mav therefore be applied to the MG equations in order 
to obtain the global smoothness of weak solutions, and hence we have proven: 

Theorem 4.1 (The MG system). Let Oq e L'^{R'^) be given. There exists a C°° 
smooth classical solution 6{t,x), of ()4.23p - (|4.27|) . evolving from 9q. 

Appendix A. Existence of weak solutions to p.ip - (|3.4p 
Here we sketch the proof of existence of global Leray-Hopf weak solutions of 
p.l|) - p.4|) evolving from ^o G L'^iR'^). We follow the general strategy used to con- 
struct weak solutions of the Navier-Stokes equations (cf. [31]). The main obstacle 
is the fact that u is obtained from 9 via a nonlocal operator of order 1. 

Denote by (— A)^/^ = A the square root of the Laplacian. Let £ C^iW^) be 
positive, with Jj^^ cj) dx ^ 1. Then (p^ = e~'^(j){x/e), for e > 0, is a standard family 
of moUifiers. We first consider the approximating system 

dtO' + {u" ■ \/)e' - AO' = -^eA^0' (A.l) 

div u" = 0, Uj = d^T^jO" (A.2) 

0^(0,.) =05, (A.3) 

where 6q — (j)^ * 9q represents the mollified initial data, and Tij are Caldcron- 
Zygmund operators. Note that ||0§||l2 < H^oIIl^ for any e > 0. 

Let s > d/2 + 1 and fix e > 0. Since A^e^ £ L^(R'^), and since eA^ gives a 
sub-critical dissipation, from standard energy arguments it follows that 

sup \\A'd'{t)\\L2 < C{e,d,cf,,T,\\9o\\L-), 
te[o,T] 

where C(e, d, </>, T, H^olU^) > is a positive constant which is finite for any T < oo. 
This a-priori estimate and a standard Galerkin approximation procedure ensures 
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the global existence of a strong H^ solution to ()A.l[) - (|A.3p . Moreover, for any e > 
we have the uniform in e energy inequality 

\r{T)fm^^) +2^ ||Vr(s)||2,(^,^ ds < ||0o|li2(R.), (A.4) 

for any T > 0, and thus 

6" is bounded in C([0, T];L^{W^)) n L^{Q, T; H^{W^)). (A.5) 

This guarantees that, up to a subsequence, 6'^ converges weakly to some function 
e L°°{0,T; L^) n L^{0,T; H^) (this convergence is weak-* in L°°{0,T;L^)). This 
does not suffice to pass to the limit in the weak formulation of (|A.ip - (jA.3p . We 
next claim that for any compact set /C C R'' we have 

dtO' is bounded in L^^^{0,T;W~^'^^{IC)). {A.6) 

Indeed, from (jA.5[) . the Gagliardo-Nirenberg inequality, and interpolation, it fol- 
lows that r is bounded in L^{0,T; L'^'^^'-'^'^^M.'^)). Since Tij are bounded from 
L2(Md) into itself, by (jXSJ) it follows that u<^ is bounded in L'^{0, T; L'^{R'^)). There- 
fore, by Holder's inequality, div(M^6'^) is bounded in 1-^/^(0, T; M^-i.2d/(2d-i)^]^d)-, 
Lastly, eA^e^ is bounded in L'^{0,T; H-^{R'^)), and A9' is a bounded family in 
^2(0, T; iJ~^(R'')). Therefore, by (jA.l[) . restricting to a compact /C, we obtain that 
dtO'^ is bounded in 

L4/3(0, T; W-^'^^ilC)) + L^{Q, T; H-^{IC)) + L^{0, T; H-\lC)), 

and hence in L'^^^{0,T; W~'^'^^^ (JC)) by the Sobolev inequality, proving (jA.6p . 

Since the injection H^{JC) into L^{IC) is compact, the injection of ^^(/C) into 
]/l/-2,2d/(2d-i)j'^-j jg continuous, it follows from the Aubin-Lions compactness lemma 
[an Theorem 3.2.1] (cf. [20]) that 

6' ^e strongly in L^{0, T; Ll^iM."^)) (A.7) 

since K. was arbitrary. Passing to the limit in the weak formulation of (|A.ip - (|A.3P 
is nontrivial only for the nonlinear term. For any (p G C^((0, oo) x M''), upon 
recalling that Uj = diTij[9], and an integration by parts in Xi, we have 

Tu" -Vip-Ou- Vip) 



'-e)u-Vip- / / d,e' T,j[9' - e] d,p - o" T,,[e' - e] d.djp 

^h + Ih + IIh. (A.8) 

Since u G LjL^, by (JA.7P and the Holder inequality it follows that It —^ as e — )■ 0. 
To obtain the convergence of 11^ and Hie ■, we claim that 

T,,[e' - 0] ^ strongly in L2(0, T; LJ^M''))- (A.9) 

The proof of (|A.9P is similar to that of (|A.7|) . Since T^ is bounded on ^^(IR'') and 
on iJi(E^), it follows from ()A.5P that 

Tij[9'] is bounded in CilO.TlL^iW^)) rM^{<d,T-H^{m.'^)) 
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Also, Tij is bounded on L2'*/(2rf-i)(M'i)^ so that we obtain Ty [u' 9"] is bounded in 
L''/3(o, T; L2<i/(2d-i)(-]jd^^^ Yhi a compact /C and a test function supported on /C. 
Applying T^j to (jA.l| . integrating against <(>, and integrating by parts, we obtain 

< CWu" e^W ^i/3 ^2d/(2d-l)\\4)\\ j^i^2,2d +C||6'^||^2^i||0||^4^y2,2d. 

In the last estimate we have also used the Holder and Poincare inequalities. The 
above proves that 

dtT.jie"] is bounded in L^'^{0,T]W~'^''^ [IC)). 

The claim (|A.9I) now follows directly from the Aubin-Lions lemma (cf. pOl I31j). 
Moreover, this shows that in (|A.9[) we have III^ -^ and //g — > as e — > 0. 

This proves that 6* is a weak solution to the limit system, i.e., p.ip - (|3.4p . By 
construction it satisfies the energy inequality, concluding the proof of existence of 
the Leray-Hopf weak solutions to p.ip - p.4p . 
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